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^ '• ABSTRACT 

O I The large H ll complex Nil in the Large Magellanic Cloud contains OB 

^ ', associations at several different stages in their life histories. We have obtained 

I ROSAT PSPC and HRI X-ray observations, Curtis Schmidt CCD images, 

\ echelle spectra in Ha and [N II] lines, and lUE interstellar absorption line 

^ \ observations of this region. The central bubble of Nil has an X-ray luminosity a 

5^ \ factor of only 3-7 brighter than predicted for an energy- conserving superbubble, 

making this the first detection of X-ray emission from a superbubble without 
a strong X-ray excess. The region Nil B contains an extremely young OB 
association analogous to the central association of the Carina nebula, apparently 
still embedded in its natal molecular cloud. We find that NllB emits diffuse 
X-ray emission, probably powered by stellar winds. Finally, we compare the 
tight cluster HD32228 in Nil to R136 in 30 Dor. The latter is a strong X-ray 
source, while the former is not detected, showing that strong X-ray emission 
from compact objects is not a universal property of such tight clusters. 



Subject headings: ISM: bubbles - Magellanic Clouds - X-ray: ISM 
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1. Introduction 

Giant H II regions contain numerous massive stars. In these regions we expect that 
supernova remnants (SNRs) and stellar winds from the most massive stars will interact 
with interstellar gas to produce shocked, hot plasma that emits X-rays. Indeed, bright 
diffuse X-ray emission has been detected in two nearby giant H ll regions: at levels of a few 
xlO^^ erg s~-^ in the Carina Nebula (Seward & Chlebowski 1982; Corcoran et al. 1994) and 
~ 10^^ erg s~^ in the 30 Doradus Nebula (Chu & Mac Low 1990, hereafter Paper I; Wang 
& Helfand 1991a). Unfortunately, these giant H ll regions are so complex that it is difficult 
to understand their observed X-ray properties. 

It has been suggested that stellar winds alone power the X-ray emitting plasma in 
the Carina nebula, since no classical SNR signatures exist, and since the diffuse X-ray 
morphology roughly follows the distribution of early-type stars (Seward & Chlebowski 1982; 
Walborn & Hesser 1982). However, interstellar absorption observations show the presence 
of high- velocity gas at AV = ±200-300 km s"^ (Walborn 1982; Walborn & Hesser 1982), 
and such large velocity offsets have otherwise only been seen in SNRs (Jenkins, Silk, & 
Wallerstein 1976; Fitzpatrick & Savage 1983). Could there be hidden SNRs producing both 
X-ray emission and high- velocity gas in the Carina Nebula? 

The diffuse X-ray emission in 30 Dor correlates well with its large shell structures, 
reminiscent of the isolated X-ray bright superbubbles described by Paper I and by Wang 
& Helfand (1991b). In these X-ray bright superbubbles, the observed X-ray luminosities 
exceed by as much as an order of magnitude the luminosities predicted by pressure-driven 
bubble models (Weaver et al. 1977; Paper I). The source of the X-ray emission has 
been proposed to be SNRs hitting the shell walls ( Paper I; Wang & Helfand 1991b) or 
mass- loaded SNRs within the bubble (Arthur & Henney 1996). We have also identified 
superbubbles that are faint in X-rays, and found that the upper limits on their X-ray 
luminosities are consistent with those expected in pressure-driven bubble models (Chu et al. 
1995, hereafter Paper HI). Can we convincingly detect X-ray emission from superbubbles 
at levels predicted by pressure-driven bubble models? 

An Einstein High Resolution Imager (HRI) observation of the core of 30 Dor showed 
two point-like sources, each having an X-ray luminosity of nearly 10^^ erg s^^ (Wang & 
Helfand 1991a). These two sources arc projected near two tight clusters of massive stars, 
R136 and R140 (Walborn & Seward 1989); however, the luminosities exceed those of stellar 
emission by at least 2-3 orders of magnitude. Using recent ROSAT HRI observations, 
Wang (1995) identified Mk34 (10" from R136) and R140 as their optical counterparts and 
concluded that these X-ray sources are Wolf-Rayet -|- black hole binaries. Do such sources 
always occur in similar tight clusters of massive stars? 

To help answer these questions, we have studied the LMC H ll complex Nil, the 
second-ranking H ll complex in the LMC (Kennicutt & Hodge 1984). Nil contains four 
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OB associations: LH9, 10, 13, and 14 (Lucke & Hodge 1970). LH9 is surrounded by a 
large superbubble, while LHIO, 13, and 14 are embedded in bright, compact H ll regions 
along the periphery of the superbubble, designated as the B, C, and E components of Nil, 
respectively (see Figure |T|). 

We have selected Nil for three reasons. First, the OB association LHIO within NllB 
contains several 03 stars similar to those in the OB associations in the Carina Nebula, but 
the surrounding H II region NllB has no high-velocity gas indicative of possible hidden 
SNRs. If the X-ray emission in the Carina Nebula comes only from hidden SNRs, we 
would expect NllB to have a much lower X-ray luminosity. Second, the large superbubble 
surrounding LH9 in Nil forms a good candidate to directly observe X-ray emission from a 
pressure-driven bubble. Third, the tight core of the association LH9, HD 32228, contains at 
least 16 early-type stars with the principal components being WC5-6 and 09.511 (Parker 
et al. 1992), analogous to R136 or R140 in 30 Dor, and so can be examined for point-like 
X-ray sources. 

Nil has recently been observed by Rosado et al. (1996) using an imaging Fabry-Perot 
interferometer to map the velocity structure in Ha and [O ill] lines. Echelle spectroscopy in 
Ha along a N-S line cutting through the main shell of Nil was performed by Meaburn et al. 
(1989). Parker et al. (1992, 1996) studied the stellar content of the central associations of 
the main shell (LH9), and of NllB (LHIO), and Walborn & Parker (1992) compared NllB 
to 30 Dor. Maps of CO tracing molecular gas in the region have been presented by Cohen 
et al. (1988), Israel & de Graauw (1991), and Caldwell (1996). 

In the next section, we describe observations of Nil in the visible, ultraviolet, and 
X-ray wavelengths that we have obtained in order to better answer the questions raised 
above and understand the dynamics of this region. In §3 we define the X-ray source regions 
and describe the analysis of the X-ray data. In §4 we derive the X-ray luminosity predicted 
by the pressure-driven bubble model. In §5 we attempt to interpret the observations of the 
various parts of the Nil complex, and finally in §6 summarize our conclusions. 

2. Observations 

We have observed Nil using a number of different techniques. The distribution of the 
denser, photoionized gas in shell walls and molecular cloud surfaces can be seen in optical 
emission-line images that we obtained using a CCD camera on the Curtis Schmidt telescope 
at the Cerro Tololo Inter- American Observatory (CTIO). In order to trace the hot gas 
produced by shocks in stellar wind bubbles or SNRs, we used pointed X-ray observations 
with the ROSAT Position Sensitive Proportional Counter (PSPC) and High Resolution 
Imager (HRI). The kinematics of the photoionized gas are revealed by high-dispersion, 
long-slit spectra obtained with an echelle spectrograph on the CTIO 4m telescope, and 
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high-dispersion UV spectra obtained with the International Ultraviolet Explorer (lUE). In 
this section, we describe these observations. 

2.1. Curtis Schmidt CCD images 

Multiple optical CCD images of two fields of the Nil complex were obtained with 
the Curtis Schmidt telescope at CTIO on 1995 January 28. A front-illuminated Thomson 
1024x1024 CCD was used; its 19 fim pixel size corresponded to 1'.'835 pixel"^. The 
field-of-view was 31'.3. Narrow-band filters of Ha (A^ = 6564 A, AA = 20 A), [S ll] 
(Ac = 6724 A, AA = 50 A), and [O ill] (A^ = 5010 A, AA = 50 A) were used to isolate 
diagnostic emission lines. Broader red and green filters (Ac = 6840 A, AA = 95 A and 
Ac = 5133 A, AA = 92 A, respectively) were used to obtain emission-line-free continuum 
images, which provide templates for continuum subtraction from the emission-line images. 
The data were reduced using the IRAFQ software package for bias subtraction, fiat-fielding, 
and sky subtraction. The images of the two fields were mosaicked after multiple frames in 
each field were shifted and combined. The total exposure times for each position in the 
final mosaicked images are 2700 s in Ha, 4800 s in [S ll], 3600 s in [O ill], 1800 s in the red 
continuum, and 1500 s in the green continuum. Figure ^(a) presents the green continuum 
image to show the location of stars, and Figure ^(b) the continuum-subtracted Ha image 
to show the distribution of ionized gas. 

The fiux calibration of the CCD images was made using observations of photometric 
standard stars. The sky-subtracted Ha fiuxes of Nil determined from the CCD images 
agree within 10% uncertainty with those determined from the photoelectrically calibrated 
PDS scans of Kennicutt & Hodge's (1986) photographic plates, kindly provided to us by 
Dr. R. C. Kennicutt. The Ha fiuxes are converted to emission measures, assuming an 
electron temperature of 10"^ K. The faintest nebulosities detected by the deep CCD images 
have emission measures of ~ 10 cm~^ pc; the brightest parts have emission measures over 
10^ cm~^ pc. 

2.2. ROSAT X-ray Observations 

We used the ROSAT PSPC and HRI to obtain X-ray observations of Nil. The PSPC 
is sensitive to the energy range of 0.1-2.4 keV, while the HRI is sensitive to 0.1-2.0 keV. 
The PSPC has an energy resolution of ~45% and an on-axis angular resolution of ~30" at 
1 keV. The HRI has a very poor energy resolution but a higher angular resolution, ~5". 
Detailed information about the performance and characteristics of these instruments can be 
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found in Pfeffermann et al. (1986) or the RO SAT Mission Description (1991). 

The PSPC observation of Nil was numbered rp900320. A total integration time of 
29.65 ks was obtained, with 16.7 ks between 1992 November 16 and 1992 December 4 
and 12.9 ks during 1993 February 14-15. The PSPC was operated in the low-gain state 
throughout the observation. The HRI observation of Nil was numbered rh900321. A total 
integration time of 27.1 ks was obtained between 1992 October 25 and 1993 February 18. 

The PSPC observation is best suited for revealing diffuse X-ray emission. To maximize 
the S/N ratio without significantly degrading the spatial resolution, we have binned the 
PSPC image, which is heavily oversampled compared to its instrumental resolution, by a 
factor of 10 to obtain a pixel size of 5" pixel"^, and smoothed the image with a gaussian 
of 0" = 5 pixels, or 25". The smoothed PSPC image in the 0.1-2.4 keV band is shown in 
Figure |^. To show the relationship between Ha and X-ray emission. Figure |^ presents an 
Ha image overlaid by X-ray contours derived by binning the PSPC image by a factor of 8 
and smoothing it with a gaussian of o" = 2 pixels, or 8". 

The HRI observation is better suited for revealing point sources. For a 6"-diameter 
detection cell, no point sources are detected in Nil and its vicinity at > 3cr levels, or 
2.3x10"^ counts s~^. This count rate corresponds to a luminosity of 4x10^^ erg s~^, using 
the conversion factor for a 5 x 10^ K plasma and an absorption column density of log = 
21.5 in Figure 10.3 of the ROSAT Mission Description (1991). 

The SNR NllL (the bright source on the right) is the only obvious X-ray source in 
the HRI image. We have binned the HRI image by a factor of 10, resulting in a pixel size 
of 5" pixel"^, and smoothed it with gaussians of a = 50" - 100". Regions with bright 
diffuse X-ray emission detected by the PSPC are also detected in the heavily smoothed HRI 
images. Since these heavily smoothed HRI images have effective resolutions comparable to 
or lower than the PSPC's resolution, and since the smoothed HRI image has lower S/N 
ratios than the PSPC image, the HRI image is not shown in this paper. Analysis of the 
HRI image of NllL will be reported in another paper. 

2.3. Long-Slit Echelle Spectroscopy of Emission Lines 

We obtained high-dispersion spectra of Nil using the echelle spectrograph with the 
long- focus red camera on the 4m telescope at CTIO. A 79 1/mm echelle grating was used. 
We used the instrument as a single-order, long-slit spectrograph by inserting a postslit 
interference filter and replacing the cross-dispersing grating with a flat mirror, so that a 
single echelle order around Ha and [N ll] AA6548, 6583 could be imaged over the entire 
slit. The data were recorded with a Tek 2048x2048 CCD using Arcon 3.6 controller. The 
pixel size was 24 /xm pixel"^, which provided a sampling of 0.082 A (3.75 km s~^) along 
the dispersion and 0"267 on the sky. Spectral coverage was effectively limited by the 
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bandwidth of the interference filter to 125 A. Coverage along the slit was 245", including an 
unvignetted field of ~ 3f5. The instrumental profile, as measured from Th-Ar calibration 
lamp lines, was 16.1 ± 0.8 km s~^ FWHM at the Ha line. This range includes variations in 
the mean focus and variations in the focus over a single spectrum. Spatial resolution was 
determined by the seeing, which was ~ 1". The exposure times were 10-20 minutes. 

The data were bias subtracted at the telescope, and later processed using the IRAF 
software for fiat-fielding, wavelength calibration, distortion correction, and removal of 
cosmic rays. Wavelength calibration and distortion correction were performed using Th-Ar 
lamp exposures taken in the beginning of each night. To eliminate velocity errors due to 
flexures in the spectrograph, the absolute wavelength scale in each spectrum was referenced 
to the geo coronal Ho; line at zero observed velocity. 

Three slit positions were observed in Nil. One E-W oriented slit was centered on the 
star Sk-66°31 (RA 04^56°'40^ Dec -66°33'33", J2000) near the center of the main shell 
in Nil. The other two slit positions were centered on the star PGMW 3021 (Parker et 
al. 1992) in NllB ~ 244"N of Sk-66°31, with one slit E-W oriented and the other N-S 
oriented. The slit positions are marked in Figure |^(b) and the echellograms are presented 
in Figure ^. The spectral features detected in the echellograms (Figure ^ include the Ha 
line and the [N ll] AA6548, 6583 lines, in addition to the faint telluric airglow features of 
OH. The [N ll] lines are too weak to be of much use. 

2.4. High-dispersion lUE Observations 

High-dispersion lUE spectra are available for two stars in Nil: SWP 15065 for 
HD 32228 (Sk-68°28) and SWP 49304 for Sk-66°31. The spectrum of HD 32228, retrieved 
from the lUE archives and presented by Chu et al. (1994), shows C iv and possibly Si iv 
interstellar absorption at ~ 270 km s~^, and S ll and Si ll at ~ 290 km s~^. The spectrum 
of Sk— 66°31, an 091 star, was obtained by us in the 16th Episode of lUE. Unfortunately, 
this spectrum is too noisy to be very useful, although by much smoothing one can find 
features that could be Magellanic C iv and possibly Si iv absorption. 

3. Analysis of X-ray Data 

To analyze the PSPC X-ray data, we first define source regions (detailed in §3.1), 
then extract the count rate and spectral energy distribution for each region. The 
observed spectral energy distribution is a convolution of the intrinsic spectrum, foreground 
absorption, and the instrumental response. The instrumental response of the ROSAT PSPC 
is well calibrated. It is possible to fold the instrumental response into plasma emission 
models (e.g., Raymond & Smith 1977) and fit the results to the observed spectral energy 



- 7- 



distribution in order to determine the plasma temperature kT, foreground absorption A^h, 
and X-ray luminosity L^. We can express the foreground absorption in terms of hydrogen 
column density A'h by using the effective cross-sections per hydrogen atom of Morrison 
& McCammon (1983). The goodness of our spectral fits is shown by gi'id plots in the 
kT — A^H plane. Using a 30% cosmic abundance that is appropriate for the LMC, we have 
made spectral fits to the observed spectral energy distribution (Figure |^) and generated 
grid plots (Figure ||). The best-fit kT, log A^h, and are listed in the top section of Table 
1. 

The number of counts received from our sources is usually too low to constrain the 
spectral fits well. Therefore, we use other independent estimates of foreground absorption 
to confine the plasma temperature fits. The various methods we use to estimate absorption 
columns are described in detail below in § 3.2. The plasma temperature kT and the 
absorption column A^h are then used to determine the intrinsic X-ray fluxes and luminosities. 
The best independently estimated log A^h, the corresponding best-fit kT, and the resultant 
Lx are listed in the bottom section of Table 1. 

When absorbing column densities are of order 10^^ cm~^ or higher, emission below 
about 0.5 keV is almost completely absorbed, even though RO SAT is sensitive down to 0.1 
keV. The intrinsic luminosity in the 0.1-2.4 keV band is then dominated by the absorption 
correction in the 0.1-0.5 keV range, but virtually no photons are detected in this range, 
giving very uncertain results. We do give the luminosities derived for both 0.1-2.4 keV and 
0.5-2.4 keV bands in Table 1. Although we used the 0.1-2.4 keV band to find conservative 
upper limits for X-ray dim superbubbles in Paper III, the high luminosities given in Table 
1 from considering this full band in regions with high absorption appear unphysical. The 
luminosities derived from the 0.5-2.4 keV band are more physically meaningful. 

3.1. Selection of Source Regions 

There is excellent correspondence between the X-ray emission and interstellar structures 
revealed in the Ha image. Frequently, enhanced X-ray emission is seen within large shells. 
Guided by the Ha image, we divide the X-ray emission from Nil into seven source regions, 
as shown in Figure ^ and listed in Table Four of these regions are associated with large 
shells. We define "Shell 1" as the main shell of Nil, encompassing LH 9. "Shell 2" is the 
faint shell structure 10' north of the main shell. "Shell 3" is the small shell structure NW 
of Shell 1 and SW of Shell 2. "Shell 4" is the shell located 10' south of the main shell. Two 
source regions are selected to include the B and C components of Nil, respectively. We 
notice faint Ha filaments correlated with diffuse X-ray emission extending NW of Shell 2 
suggesting a blowout into a lower density region which we define as the "NW diffuse" source 
region. The X-ray data of these regions are analyzed individually. 
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3.2. Independent Estimates of Absorption Columns 

Soft X-rays are absorbed by all phases of the interstellar medium (ISM). We have used 
two methods to estimate the ranges of absorption column densities. As described below, 
both methods suffer from large uncertainties. 

In the first method we limit the possible range in absorption column A^h by using 
previous observations of atomic and molecular gases and our current observations of ionized 
gas toward Nil. An H I survey of the LMC by Rohlfs et al. (1988) has shown a range of 
LMC H I column density of (1.8-2.9) xlO^^ cm~^ in the four pointings near Nil, with a 
beam size of 15'. A Galactic foreground H I column density of Ah = (5.6 ± 2.0) x lO^'^ cm~^ 
was derived by Schwering & Israel (1991) using the large-scale H I survey by Cleary, Heiles 
& Haslam (1979). 

A high-resolution ESO-SEST survey shows a molecular cloud complex roughly encircling 
the Nil complex. Except for NllB, X-ray emission is not observed from any regions 
containing CO, possibly due to obscuration by the molecular clouds. The clouds show a 
peak CO intensity, Ico = 20 K km s~^ (Israel et al. 1993). To translate this into a column 
density of Il2, we scale the galactic CO/H2 conversion ratio Amw = ^n-z/^co = 2.3 x 10^° 
cm~^ K^^ km^^ s (Bloemen 1989) by the LMC metallicity of 1/3. The observed CO 
intensity then correspond to a of ~ 1 x 10^^ cm~^. 

Using our calibrated CCD Ha images, we find emission measures of 4 x 10^ cm~^ pc 
in NllB and less than 1 x 10^ cm~^ pc in shells and other low-surface brightness areas. 
Assuming the depth of the Ha emission region is as large as the width of the emission 
region, we estimate ionized hydrogen column densities of ~ 3.5 x 10^^ cm~^ in NllB and 
~ 2 X 10^° cm~^ in shells and low-surface brightness areas. 

By combining these values we can estimate upper and lower limits to the absorbing 
hydrogen column density. For the lower limit, we assume that Nil is in front of all 
the LMC H I and that no more than half of the ionized gas in Nil contributes to 
the absorption column, so that the lower limit is the sum Galactic Nm + | LMC Ann. 
We use the appropriate H ll column densities for the different regions as given above. 
For NllB this yields a lower limit of log Ah = 21.3, and log Ah = 20.9 for the other 
regions. For the upper hmit, we examine the CO image of Caldwell (1996). In NllB, 
the only X-ray emitting region that coincides with strong CO emission, we use the sum 
Galactic Ahi + |( LMC Ann + LMC Am + 2 LMC AhJ as an upper estimate, while in 
regions lacking CO emission we neglect the H2 contribution entirely. This gives an upper 
limit of log Ah = 21.6 for NllB, and log Ah = 21.3 for the other regions. 

The second method to estimate absorption column uses optical extinction of stars 
embedded in the X-ray source regions. An overall average reddening of E{B — V) — 0.05 
mag has been found in the stars of LH 9, the OB association in the main Nil shell, and 
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E{B - y) = 0.17 mag in LH 10 in Nil B (Parker et al. 1992). Brunet (1975) has reported 

an average Galactic foreground reddening of E{B — V) = 0.07 mag. However, small scale 
variations in Galactic reddening ranging from 0.0 to 0.15 mag toward the LMG have 
been found by Oestreicher et al. (1995) using a more complete sample of 1500 Galactic 
foreground stars. We adopt a Galactic gas-to-dust ratio of Nii/E{B — V) = 4.8 x 10^^ 
cm~^ mag~^ (Bohlin et al. 1978), and an LMG gas-to-dust ratio of 2.4 x 10^^ cm~^ mag~^ 
(Fitzpatrick 1986). Adopting the reddening values of Parker et al. (1992) we estimate 
20.3 < logA^H < 21.1 for the main shell and 21.1 < logA^n < 21.6 for NllB. By the same 
arguments, we can make an independent estimate of 20.9 < logA^n < 21.6 using the LMG 
intrinsic reddening map of Oestreicher & Schmidt-Kaler (1996). 

These values agree well with our estimates from the emission-line observations of 
the ISM described above. We use a combination of the two estimates to constrain the 
absorption column densities, taking the lower of the two upper limits and the upper of the 
two lower limits to give the most likely range of A^h- The range of A^h that we have actually 
used for each region is hsted in the bottom section of Table 1. 



4. X-ray Emission from Pressure-Driven Bubble Model 

To predict the X-ray luminosity from a pressure-driven bubble we follow the procedure 
described in Paper I, as corrected in Paper III. By assuming that conductive evaporation 
determines the interior structure of the bubble (Weaver et al. 1977), we can compute the 
X-ray luminosity if we know the radius of the bubble in pc, Rpc, its expansion velocity in 
km s~^, V5, and the emission measure in Ha on a tangential path through the shell in cm~^ 
pc, EM. If we assume that the thickness of the shell, AR <^ R, the electron temperature 
Te ~ lO'^, and the mean atomic mass of the ambient medium /la = (14/ll)mij, we can 
use equations (7-10) from Paper III to derive the X-ray luminosity directly from the three 
observables, 

~ (6.7 X 10"^ erg s-^)iIEM''l' RfJ'vl'\ (1) 

where ^ is the metaUicity, which we take to be 0.3 in the LMG, and / is a slowly varying 
function of the interior temperature of the bubble which we can take to have value 7 ~ 2. 
The external density into which the bubble is expanding is then 

no ~ (22 cm.-^)EM^I^vf''R-^l\ (2) 

the age of the bubble in Myr is 

= 0.59i?pe/^^5, (3) 
and the mechanical luminosity Lmech from stellar winds and SNRs driving the bubble is 



Lmech ^ (8 X lO^'^ erg s-')EM^'\T R%\ 



(4) 
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The kinetic energy K of the shell in a pressure driven bubble can be derived by noting 
that the thermal energy of the interior is f/ = (5/ll)Linech^ (Weaver et al. 1977), and that 
K + U = 50/77Lniecht (Mac Low & McCray 1988) , so 

K = (15/77)L^echt = (2.9 x lQ^\igs)EM"\l'^ R%\ (5) 

(Alternatively, the final result can be derived by taking K = {l/2)MVg, where the mass of 
the shell M = (4/3)7ri?Va^o-) 

5. Interpretation 
5.1. Shell 1 

In Paper I we showed that some LMC superbubbles are X-ray bright, by which we 
mean that they have X-ray luminosities much higher than predicted by the pressure-driven 
bubble model, while in Paper III we placed upper limits on the X-ray luminosities from 
several other superbubbles that were close to the prediction of the model. Here we derive 
the X-ray luminosity of Shell 1 and show that it is only slightly X-ray bright, with a 
luminosity a factor of 3-7 brighter than predicted by pressure-driven bubble model. 

X-ray emission from the central region of Nil lies within the optical shell seen in Ha. 
Performing spectral fits using free parameters in this region shows a 99% confidence range 
in absorption and plasma temperature of logA^H =20.5-22.0 and kT =0.1-0.4 keV (see 
Table |l|). This absorption column density is consistent with the values we derived in § 3.2. 
If we constrain the column density to those values, we can derive an X-ray luminosity of 
(3.2-4.0) X 10^^ erg s~^ in the 0.5-2.4 keV band, assuming a distance of 50 kpc to the LMC 
(Panagia et al. 1991). 

The velocity structure of the main shell is quite complex. Figure ^(a) shows an 
echellogram from an E-W slit across the central region, while Figure ^ shows velocity 
profiles at four different points on that slit. The broadest splitting seen is some 120 
km s~^ (Fig. 0a), suggesting an expansion velocity as high as 60 km s~^. Meaburn et 
al. (1989) observed the Ha line with an echelle spectrograph along a N-S line passing 
through HD 32228, crossing our slit just west of center. They find as many as five velocity 
components in the central region, which they attribute to multiple expanding shells in a 
cellular structure. Rosado et al. (1996) performed Fabry-Perot observations of the entire 
Nil complex in both Ha and [O ill]. They also found multiple components, but concluded 
that the main shell could best be described as undergoing a general radial expansion at 
45 km s~^, with lower density interior gas undergoing more rapid, cellular expansion. 

We now can find the luminosity predicted by the pressure-driven bubble theory. The 
typical emission measures of filaments around the edge of the bubble are in the range 
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600-900 cm~^ pc. We take the expansion velocity to lie in the range 45-60 km s~^. To find 
the physical radius of the bubble, we again use a distance to the LMC of 50 kpc (Panagia 
et al. 1991), so that 1' = 14.5 pc. We then estimate the X-ray emitting bubble to have a 
radius R =55-65 pc. A pressure-driven bubble with these parameter ranges is predicted by 
equation (||) to have an X-ray luminosity = (0.6 — 1.2) x 10^^ erg s~^. This is a factor 
of 3-7 lower than the X-ray luminosity derived from the observations using constrained 
column densities: the main shell of Nil is somewhat X-ray bright. 

The excess X-ray emission is most likely produced by an intermittent source, since 
Nil represents an intermediate case between the bubbles we described in Paper III that 
show no excess, and those we described in Paper I with excesses of more than an order 
of magnitude. Arthur & Henney (1996) suggested that an SNR going off inside a bubble 
with embedded clumps could produce the excess. They ignored thermal evaporation of 
the clumps, suggesting that the evaporation could be suppressed by magnetic fields. Tao 
(1995) has shown, however, that fields cannot be tangled by dynamical means sufficiently 
to strongly suppress evaporation, so clumps would probably evaporate before being swept 
up by a SNR. Alternatively, we have suggested in Paper I that the excess X-ray emission 
is being produced as interior SNRs expanding in the low-density gas hit the ionized inner 
edge of the swept-up shell. 

Using the same estimates of the observed parameters as above, we can examine the 
physical parameters of this bubble using equations (@-§). We find that the ambient density 
is 0.1-0.3 cm~^, a reasonable interstellar value. We note that there are clearly regions 
of far higher density in the ring around the main bubble, most noticeably traced by CO 
observations showing a ring of molecular clouds roughly coincident with the regions of 
highest Ha intensity (Caldwell 1996; Israel & de Graauw 1991). However, these regions 
appear to have a low enough filling factor to not be dynamically important. 

The mechanical luminosity driving the bubble (eq. [^) lies in the range 
-^^mech = (3 — 6) X 10'^'' crg s~^. Usiug only the 18 0-type stars in LH9 that have 
been spectroscopically classified by Parker et al. (1992) and using the empirical mass 
loss rates and stellar wind terminal velocities appropriate for the spectral types (Conti & 
Underbill 1988), we derive an integrated wind mechanical luminosity of 3.3x10'^'' erg s~^. 
It is conceivable that the WC5-6 star in HD 32228 and other unclassified early type stars in 
LH9 could supply the Lmech required to drive shell 1. 

The total input energy E = L^echt in such a bubble is then (0.6-1.2) x 10^^ erg, just 
allowing the small number of supernovae suggested by Rosado et al. (1996) from detailed 
modeling of the IMF. Note that the possible hidden SNR producing the excess X-ray 
emission would not yet have fully transferred its energy to the shell, and so would currently 
make only a small contribution to the quoted input energy. The kinetic energy of the shell, 
from equation (^), is i^' = (1 — 2) x 10^*^ erg. This model is definitely inconsistent with the 
dozens of supernovae suggested by Meaburn et al. (1989) based on the overall SN rate in 
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the LMC and the size of the shell. 

The age of the bubble is more problematic. From equation and the above estimates 
of radius and expansion velocity, we find ages of 0.5-0.8 Myr, well below the age of the 
central OB association. Another way of looking at it is that this bubble is expanding 
much too fast to have been produced by such an old OB association. Oey (1996) has 
demonstrated that this same phenomenon occurs in a number of superbubbles where the 
age of the central OB association has been measured. A possible solution to this problem 
may be that the expansion of the superbubble in a relatively low density medium only 
occurs after it has broken out of its parent molecular cloud, well after the formation of the 
central OB association. The association LH 10 in the nebula Nil B, which we discuss next, 
may in fact be an example of just such an association that has formed, but has not yet 
broken out into the low density medium and produced a clear supershell. This suggestion 
supposes that the boundary between dense molecular cloud and diffuse ISM is significantly 
steeper than the power law that Oey (1996) examined. Further modeling of the process 
of superbubble breakout from molecular clouds appears useful. 



5.2. NllB 

An unconstrained spectral fit to this region gives a range in column density 
logA^H = 20.3 — 22.0 and temperature kT = 0.1 — 1.0 keV. Constraining the absorbing 
column density independently, we find more reasonable values of logA^H =21.3-21.6 and 
kT =0.4-0.5 keV. This corresponds to an X-ray luminosity = (0.3 — 0.8) x 10'^^ erg s~^. 

In order to estimate the X-ray emission contributed by OB stars in the central cluster 
LH 9, we first convert V band fiuxes of the 24 O and B stars located in a 4 arcmin^ region 
centered on HD 32228 into bolometric magnitude (Mboi) by applying reddening (Parker et 
al. 1992) and bolometric (Chlebowski & Garmany 1991) corrections. We then convert the 
summed Lboi to X-ray luminosity by the relation = (1.4 ± 0.3) x 10~^Lboi, which is good 
for stars of types between 03 and A5 (Pallavincini et al. 1981; Chlebowski, Harnden, & 
Sciortino 1989). The overall stellar X-ray luminosity is 5x10^^ erg s~^, which is 6-17% 
of the observed X-ray luminosity in Nil B. 

Figures ^(b) and^(c) show two echellograms of the Nil B region. The two slit positions 
are in a cross roughly centered on Nil B, as shown in Figure ^(c) and described in § 



They show well-resolved, remarkably Gaussian line profiles, with FWHM of 30-40 km s~^, 
far broader than the instrumental profile. Figure ^ shows an example of one of these profiles 
at one of the few places where a high-velocity shoulder might be visible on the line. Other 
than this shoulder, the lines appear quite featureless, and completely lack high-velocity 
emission at velocities greater than 100 km s~^. 



The Carina Nebula resembles NllB in stellar content, both containing several 03 
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stars. A major difference between them is that the Carina nebula shows an expanding shell 
structure in velocity and shows high-velocity gas at velocity offsets of several hundred km 
s~^ (Walborn 1982; Walborn & Hesser 1982), while NllB appears to have neither of these 
two dynamic features. The heating mechanism for the X-ray emitting gas in Carina could 
be either stellar winds or SNRs. 

The lack of a visible ring structure in NllB suggests that it contains no large cavity of 
hot gas. This makes it difficult to hide an SNR in the region. A supernova going off in a 
region with density greater than 0.1 cm~^ will generate classical SNR signatures (Tomisaka, 
Habe, & Ikeuchi 1981). We can find a lower bound to the density in NllB by assuming 
that it is a uniform, ionized sphere. The emission measure in Ha is ~ 4 x 10'^ cm^® pc, and 
the diameter of the region is 27 pc, yielding an electron density of Ue — 38 cm~^. An SNR 
expanding into gas this dense would produce fireworks, in the form of the classical SNR 
signatures, none of which are observed. LH 10 contains between three and six 03 stars 
(Parker et al. 1992), suggesting that the cluster is extremely young, and adding further 
weight to the suggestion that no SNR has occurred yet. 

If stellar winds do produce the observed X-ray emission, however, the question arises 
of why no clear stellar wind bubbles can be seen in NllB. A pressure-driven bubble would 
have a radius of (Weaver et al. 1977) 

R=i28pc)iM,vl,,n,HI)'/', (6) 

where Mg is the mass-loss rate of the stellar wind in units of 10^^Mq/jt, ^2000 is the stellar 
wind velocity in units of 2000 km s~^, and is the age of the bubble in units of 10^ yr. 
Substituting a density of Uq = 38 cm~'^ and setting the other values to unity gives a bubble 
of radius 14 pc, or nearly an arcminute at the distance of the LMC. Indeed, Rosado et al. 
(1996) suggest a structure of overlapping bubbles based on their Fabry- Perot observations 
of the region in Ha and [O ill] . They found moderately broad Ha lines, consistent with our 
echelle observations, as might be expected from many overlapping bubbles, and found that 
the narrower [O III] lines showed splitting indicative of small bubbles in some regions. Our 
Ha image also shows some structure suggestive of overlapping bubbles; this can be seen 
even more clearly in the Ha image of Nil B shown in Walborn & Parker (1992). We have 
been awarded time to perform high-resolution imaging with the Hubble Space Telescope 
that may resolve this structure more clearly. AXAF imaging may reveal whether X-ray 
emission is confined to the interiors of such small wind-blown bubbles. 



5.3. Other Shells 

Shell 2 lies on the north edge of Nil B. In Ha emission, the shell appears elongated 
in the N-S direction, with a diameter of approximately 5' in that direction. However, 
the combination of the X-ray image and echelle spectroscopy by Meaburn et al. (1989) 
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shows a different picture. X-rays are only emitted in the southern third of the bubble (see 
Figure |^). Spectral fits with constrained absorption give a kT range 0.3-0.4 keV and an 
X-ray luminosity of (0.95 — 1.1) x 10^^ erg s~^ in the 0.5-2.4 keV band. This region also 
turns out to show a line splitting of ~ 20 km s~^ in the echelle spectrum. The region 
without X-ray emission also shows no line-splitting in Ha. Thus the X-ray emission is 
associated with a kinematically distinct structure of unknown origin. 

Shell 3 lies adjacent to Shell 2, north of the main bubble and NW of Nil B. A kT range 
0.1-0.5 keV and an X-ray luminosity range of (0.7-0.9) xlO^^ erg s~^ are derived from the 
constrained iVn fits. There is no clear boundary between the X-ray emission in the southern 
region of Shell 2 and that of Shell 3, although the Ha morphology suggests that they are 
separate structures. The X-ray emission from Shell 3 actually extends further to the west, 
possibly even merging with the plume of emission north of the SNR Nil L. Spectroscopy of 
the region would be needed to further understand the relation among these structures. 

Shell 4 lies south of the main bubble. Even the unconstrained fit shown in Figure ^ is 
quite close to the independently derived constraints on A^h- Spectral fits with constrained 
absorption give a kT range 0.3-0.5 keV and an X-ray luminosity range of (0.97-1.2) x 10^^ 
erg s~^. It is perhaps most remarkable in that it appears quite similar in Ha morphology to 
the shell immediately to its E, and yet has detectable X-ray emission, while the other shell 
does not. In the absence of any kinematical information from spectroscopy, we can only 
suggest that intermittent events such as hidden supernovae might produce such a pattern. 



5.4. HD 32228 

X-ray emission is observed in the PSPC map peaking at the position of the bright 
cluster HD 32228 (RA 04'^56""34?5; Dec -66°28'25"; J2000) in the center of the main Nil 
shell (marked in Figure 1(a)). However, no X-ray emission is seen in the HRI observation 
over 3 0" or 2.3 x 10~^ counts s~^. This indicates that the X-ray luminosity from any point 
source associated with HD 32228 is below an upper limit of Lx ~ 4 xlO^^ erg s~^. 

The observed PSPC count rate within the 4 arcmin^ region is 9.50 x 10"^ counts 
s~^. If we take the absorbing column density and local plasma temperature to be roughly 
the same as elsewhere in Shell 1 [kT ~ 3 x 10^ K, and log Ah ~ 20.5, from Table 1), the 
observed X-ray luminosity is = 4.8 x 10'^^ erg s~^. Following the same procedure used 
in § 5^, we estimate the stellar X-ray luminosity to be ~ (1.0 ± 0.5) x 10^^ erg s~^. 
Therefore, the stellar X-ray emission accounts for only ~ 30% of the observed L^, and the 
remaining emission is from the hot gas. HD 32228 does not appear to contain a compact 
object similar to those observed near R136 and R140 in 30 Dor (Wang 1995), showing that 
such X-ray sources are not universally associated with tight clusters such as this. 
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5.5. NW region 

The Ha image shows faint filaments extending above Shell 2 to the NW, associated 
with diffuse X-ray emission, as well as with an unidentified X-ray point source, all lying 
within the box labeled NW diffuse in Figure ^. The spectrum of the point source suggests 
it is a Galactic foreground star. Its point nature has also been verified by its spatial profile 
in the HRI image. Therefore, X-ray emission from the point source is subtracted from the 
total PSPC flux. Spectral fits to the remaining emission with constrained absorption give a 
kT of 0.3 keV in this region, and an X-ray luminosity range of (0.5 — 0.6) x 10^^ erg s~^. 
Similar diffuse X-ray emission has been observed outside the major structures visible in Ha 
in N44 (Paper II) and 30 Dor (Chu 1996). 

5.6. Nil C 

The X-ray emission in this region is dominated by a weak, point-like source. Spectral 
fitting suggests a temperature of around 1 keV and a low absorption column density of 
logA^H = 20.4. These properties are consistent with a galactic stellar source in the solar 
neighborhood. Unfortunately, we cannot confirm the point nature of this source, as it has 
luminosity below the detection limit of our HRI observation. 

6. Conclusions 

1. The main bubble of Nil, surrounding the OB association LH 9, has an X-ray 
luminosity of (3 — 4) x 10^^ erg s~^, a factor of 3-7 times higher than the prediction of an 
energy-conserving, pressure-driven superbubble model. This could be the result of an SNR 
hitting the inner edge of the shell (Paper I), or of a mass- loaded SNR in the interior (Arthur 
& Henney 1996). Both of these processes should produce intermittent and time-dependent 
excesses of X-ray emission over the prediction of the pressure-driven model. Indeed, the 
X-ray emission from the main bubble of Nil is only slightly above the predicted value, 
and the bubbles of Paper HI all have upper limits on their luminosities consistent with the 
prediction. It remains true that no superbubble has yet been positively detected emitting 
with luminosity within a factor of three of the prediction of the pressure-driven model, 
although the limits are narrowing. 

2. Diffuse X-ray emission is detected from Nil B, a region analogous to the Carina 
nebula with multiple 03 and other massive stars. The luminosity greatly exceeds that 
expected from the stars themselves, as in Carina. Although no single bubble is observed, 
echelle observations show strongly supersonic motions. Ha images show hints of overlapping 
shell structures, and the Fabry-Perot observations of [O ill] by Rosado et al. (1996) show 
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line splitting in some regions. No gas moving at velocities over 60 km s~^ is seen however, 
suggesting a lack of SNRs. Thus we conclude that the diffuse X-ray emission could be 
powered solely by stellar winds from the central massive stars of LH 10. 

The application to Carina is, unfortunately, less than incisive. An absence of emission 
from NllB in excess of stellar emission would have been strong support for the suggestion 
that hidden SNRs produce the X-ray emission in Carina, because the clusters are otherwise 
similar. However, the presence of bright X-ray emission from NllB docs not argue against 
hidden SNRs in Carina, because of the presence of high-velocity gas in Carina and the 
absence of such gas in NllB. 

3. The dynamical age of the main shell of Nil is under 1 Myr, significantly shorter than 
the age derived by Parker et al. (1992) for the central cluster, LH 9. This same problem 
has been found for other superbubbles by Oey (1996). We propose that the solution to 
this problem is that during the initial evolution of OB associations, they remain cocooned 
within their high- density parent molecular clouds, and do not begin generating the observed 
superbubbles until after they have broken out of the dense clouds into the surrounding 
diffuse ISM. Nil B may be an example of a young OB association with energetic winds that 
has not yet broken out of a high-density region, since despite the presence of 03 stars and 
X-ray emission, very little high-velocity gas is observed, while CO is still observed towards 
this region. 

4. The tight cluster HD 32228 in LH 9 is not a strong X-ray source, despite its 
resemblance in stellar density and size to the clusters R136 and R140 in 30 Dor. This 

suggests that the strong X-ray emission observed from R136 and R140 is not an intrinsic 
property of tight stellar clusters, but is rather generated by a compact object that has 
formed there, as suggested by Wang (1995). 

We thank G. Garcia-Segura for participation in early stages of this project. SDP 
acknowledges the hospitality of the Max-Planck-Institut fiir Astronomic, and support from 
a National Science Foundation Graduate Fellowship. This research was partially supported 
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Figure Captions 

Fig. 1. — (a) Green continuum image of Nil from the Curtis Schmidt, showing the 
distribution of stars. The OB associations and HD 32228 are marked, (b) Ha Schmidt 
image of Nil to show the distribution of ionized gas. The OB association LH 9 hes in 
the center of the main shell; LH 10, LH 13, and LH 14 lie in Nil B, Nil C, and Nil E, 
respectively. The compact cluster HD 32228 lies just east of the center of the main shell, at 
RA 04^56"'34?5, Dec -66°28'25" (J2000). The three echelle slit positions are marked. 

Fig. 2. — X-ray image of the Nil region. The PSPC image was binned by a factor of 10 
and smoothed with a gaussian of cr = 25". The RA and Dec are in J2000. The contours are 
at 10%, 15%, 20%, 25%, and 30% of the peak brightness of the supernova remnant NllL at 
04^54"^48?5, -66°25'45". 

Fig. 3. — Contours of X-ray emission superimposed on our Curtis Schmidt Ha image. This 
X-ray image was derived by binning the PSPC image by a factor of 8 and smoothed by a 
gaussian of a = 8". The contours are 10%, 13%, 16%, and 20% of the peak brightness of the 
supernova remnant NllL at 04^54™48?5, -66°25'45". The RA and DEC are in J2000. 

Fig. 4. — Echellograms of Nil. (a) E-W oriented slit centered on Sk -66°31 close to the center 
of the main shell, (b) E-W oriented slit 244" N, in Nil B, centered on the star PGMW 3021. 
(c) N-S oriented slit also crossing the star PGMW 3021 in Nil B. The profiles shown in 
Figures |^ and |^ are drawn from the regions bounded by the thin black lines. Unvignetted 
slit lengths are 3^5. 

Fig. 5. — Observed spectral energy distributions and the best model fits without constraints 
on A^H- 

Fig. 6. — gi'id plots for the spectral fits shown in Figure ^. The three contours correspond 
to 68%, 90%, and 99% confidence levels. Temperatures are given in keV. 

Fig. 7. — Velocity profiles at the four places designated by the thin black lines in Figure §(a), 
moving from left to right. 

Fig. 8. — Velocity profile at the place designated by the thin black lines in Figure §(c). 
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Table 1. Observed X-ray properties of Nll'^ 



Regions Shell 1 Shell 2 NW diffuse Shell 3 Shell 4 Nil B 
Free N^^ fits 

kT (keV) 0.1-0.4 0.1-0.5 0.1-0.6 0.1-0.5 0.3-0.5 0.1-1.0 

logA^H 20.5-22.0 20.4-22.0 20.0-22.0 20.2-22.0 20.0-20.8 20.3-22.0 

logL^ (0.1-2.4)^ 35.57-39.36 35.00-38.85 34.61-38.52 34.78-38.65 34.92-35.28 34.39-38.26 

logL^ (0.5-2.4)^ 35.39-38.11 34.85-37.60 34.46-37.27 34.62-37.40 34.77-34.98 34.29-37.01 

Fixed A^H fits 

kT (keV) 0.31-0.33 0.34-0.36 0.30-0.32 0.27-0.29 0.26-0.27 0.36-0.51 

logA^H 20.9-21.1 20.9-21.1 20.9-21.1 20.9-21.1 20.9-21.1 21.3-21.6 

logL^ (0.1-2.4)'^ 35.73-35.84 35.18-35.28 34.90-35.02 35.09-35.22 35.27-35.38 34.69-35.09 

logL^ (0.5-2.4)=^ 35.51-35.60 34.98-35.07 34.67-34.77 34.83-34.94 34.99-35.09 34.53-34.89 



^In units of erg s ^, with the energy band given in keV. 
''Parameter ranges give 99% confidence levels. 
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